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Abstract 



I review the effects of dust, both in absorption and in emission, on the spectral 
energy distributions (SEDs) of Local star-forming galaxies. The energy balance 
between the stellar light absorbed by dust at UV-optical-nearlR wavelengths and the 
energy re-emitted by dust in the far infrared (FIR) shows that the amount of dust 
O ■ extinction affecting the galaxy can be predicted with simple prescriptions applied 

I to the UV-optical data. The implications of these results for high redshift galaxies 

^ ■ are discussed. Arguments can be given to support the view that Local starbursts 

, are representative of high-redshift (z>2), UV-detected, star- forming galaxies. If this 

is the case, the high redshift FIR emission will be generally undetected in sub-mm 



■ surveys, unless (1) the bolometric luminosity of the high-z galaxies is comparable 



to or larger than that of Local ultraluminous FIR galaxies and (2) their FIR SED 



P ' contains a cool dust component. 
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1 Introduction 



A global measure of the impact of dust obscuration on the light emerging from 
galaxies is provided by the cosmic background light spectrum. The background 
light above «i40 /im ( [Hauser et al., 19^ ; [Fixsen et al., 199"8| ) contains between 



50% and >70% of the bolometric emission from galaxies integrated over all 
redshifts (Figure 1). This wavelength range is where most of the stellar energy 
absorbed by dust is re-emitted. 

While the redshift-integrated effect of dust on galaxy light is somewhat con- 
strained, the relative importance of dust obscuration as a function of redshift 
is not well established, with the possible exception of the Local Universe. In 
nearby galaxies, the 8-120 /xm FIR emission measured by IRAS represents only 
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Fig. 1. The background light emission from the ultraviolet to the far infrared, given 
in terms of energy as a function of wavelenght. The UV data point at 0.2 /im is 



from Armand, Milliard &: Deharveng (1994 ), the optical data are from the Hubble 



Deep Field, and the K-band data point from ground-based surveys (Pozzetti et al. 
1998). The lower limits in the FIR region are from IRAS in a no-evolution scenario. 



The upper limits and the triangles at A >60 /xm are from COBE-DIRBE ([Hausei] 
et al., 1998). The smooth curves are the analytical representation (dot-dashed line) 
and the 1 a envelope (dotted lines) of the background detected by COBE-FIRAS 
([Fixsen et al., 199^ ). Adapted from Pei, Fah fc Hauser (1999] ). 



~25% of the bolometric output from galaxies ( [Soiter fc Meugebauer, 19^ . 
However, when compared with the galaxy UV output, the energy in the IRAS 
bands represents between 1/2 and 2/3 of the total UV+FIR emission in Local 
galaxies. Regions of active star formation emit the bulk of their energy in the 
UV, owing to the presence of young, massive stars, and the dust-reprocessed 
UV light is the major contributor to the FIR emission in the IRAS wavebands. 
Thus, between half and two-third of the star formation in the Local Universe 
is obscured by dust. 

Observational evidence points to increased dust obscuration in more distant 
galaxies ( (Jalzetti fc Heckman, 195^ ; [Steidel et al., 195^ . Part of this effect is 
due to the fact that standard observing techniques access bluer and bluer, and 
thus more and more dust- sensitive, restframe wavelength regions for increasing 
redshift. Using dust-correction methods tested in the Local Universe ( |Calzett 



"199^ ; IMeurer, Heckman fc Calzetti, 1999| ; [Calzetti et al., 2000|) , the UV light 
from high redshift star-forming galaxies has been shown to be possibly ob- 
scured at the 60%-80% level. If confirmed by further evidence, this result has 
a major impact on a number of properties of the high-z Universe, including 
the star formation rate (SFR) estimates ( |Steidel et al., 1999| ). Moreover, the 
selective absorption of dust as a function of wavelength is likely to be the main 
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Fig. 2. The average star formation rate density of the Universe in Mq yr~^ Mpc~^ 
as a function of redshift, before (left panel) and after (rigth panel) correction for 
the effects of dust obscuration. Circles indicate restframe UV data at 0.28 (for 
z<2, [Lilly et al. (T996| ); [Connolly et al. ( 19971) ) and at 0.17 (for z>2, [Steidel et] 
al. (1999|)). Triangles mark the position of the restfra me Ha measurements (palleg^ 
et al., 19951 ; ICronwah, 199^ ; [Tresse k Maddox, 1998| ; [Glazebrook et al., 1999| ; [Ya: 
1999|) . Squares mark the position of mid/far-IR data, from ISO (z<l, Florei 



et al., 



et al. (1999 | ) ) and from SCUBA ([Hughes et al., 1998| ; [Barger, Cowie k Richards 



2000). 1 a error bars are reported for all the data. 



reason for the discrepant results given by different SFR indicators. Figure 2 
shows the SFR density of the Universe as a function of redshift ( pVIadau et al 



1996| ) measured from restframe UV, Ha, and FIR emission, before and after 
dust correction. Only after the effect of dust is factored in the estimates, the 
three different SFR indicators are brought back into agreement. 

Quantifying the impact of dust on the light emitted by the stellar component 
of galaxies bears not only on our understanding of the evolution of the star 
formation, but also on a host of other aspects of galaxy evolution, including 
the relationship between star formation and AGN. Various authors ( Poyle & 



'I'erlevich, 1998 ; [Shaver et al., 19TJ^) have pointed out the close resemblance 
between the SFR density of galaxies and the space density of quasars (the 
latter measured in the radio) as a function of redshift. The comparison be- 
tween the two trends is, however, strongly dependent on the adopted dust 
corrections for the SFR of galaxies, and may or may not result in a similarity. 
For instance, corrections like those shown in the right-hand panel of Figure 2 
imply a relatively constant SFR density for z>2.5, in a regime where the space 



density of quasars is declining for increasing redshift ( [Shaver et al., 19"^ ) 
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2 Why Is Dust Elusive? 



Paraphrasing Witt, Thronson & Capuano (1992), the derivation of the opacity 
of a galaxy "depends sensitively upon the adopted [dust] geometry, largely 
due to the important effects of scattering and to the relative fraction of stars 
that are lightly obscured" . The effects of dust obscuration and reddening are 
very hard to quantify in galaxies, especially if the only data available are 
in the UV-to-nearlR, i.e. in a wavelength range where dust manifests itself 
predominantly by absorbing or scattering the stellar light. Three main factors 
conspire to make dust an elusive component of galaxies: (1) the lack of obvious 
emission or absorption features below the K band (except, in some instances, 
for the 0.22 fim absorption 'bump'); (2) the grey obscuration often produced 
by the combination of dust geometry, varying optical depth, and scattering; 
(3) the age-dust degeneracy. 

The largest majority of the methods for measuring the opacity /reddening of a 
galaxy use the light from the galaxy itself as the 'illuminating source' for the 
dust. However, the stellar populations of galaxies are complex mixtures of stars 
of all ages. A dust reddened population often resembles an older population. 
The dust-age degeneracy is an important effect when broad band colors or 
low resolution spectroscopy only are available for the object of interest and 
stellar features are not identifiable. For instance, an optical attenuation Ay=l 
changes the galaxy colors to mimic an age "increase" of a factor of ~5 in a 
stellar population younger than 100 Myr. Conversely, the colors of a relatively 
old stellar population will not be easily discriminated from the colors of a 
young but reddened population. 

This problem can be minimized in starburst galaxies, where the UV-to-nearlR 
output is dominated by the recently formed stars. In particular, massive stars 
dominate the UV emission, and their spectra are in the Rayleigh- Jeans regime 
above the Lyman discontinuity at 0.0912 /im. Thus the shape of the UV spec- 
trum of starburst galaxies is relatively constant over a fairly large range of 
ages and physical parameters ( [Leitherer fc Heckman, 1995|) ; the non- ionizing 
photons are less age-sensitive than the ionizing photons, and the nebular emis- 
sion lines from a starburst will fade before appreciable changes in the shape 
of the UV continuum can be observed. In summary, the starburst populations 
in galaxies can provide a relatively standard 'illuminating source' for studying 
the effects of dust in these objects. 
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3 Measuring Dust Obscuration in UV-Selected Starbursts 



The ionizing and non-ionizing stellar light absorbed by dust is re-emitted in the 
FIR, and measuring dust reddening and obscuration in galaxies requires multi- 
wavelength information. Figure 3 shows the SEDs of two starburst galaxies, 
covering from the UV (~0.12 /im) to the FIR (~200 fim). The two SEDs are 
of a metal-rich, FIR-luminous galaxy (NGC6090) and a low-metallicity. Blue 
Compact Galaxy (Toll924— 416). As expected, the metal-rich galaxy is also 
dusty, with about 10 times more energy emerging in the FIR than in the UV- 
optical range. In the low-metallicity galaxy, the energy share between FIR and 
UV-optical is about 50% each. Although the two SEDs are clearly different, 
the "action" of dust on the stellar light emerging from the two galaxies can 
be parameterized in a simple manner, at least in the wavelength range 0.12- 
2.2 fim ( Calzetti, Kinney fc Storchi-Bergmann, 1994 , 1996 ; [Calzetti, 1997 ; 



(Jalzetti et al., 2000| ): 

F„(A) = F,(A) 10-°-4^'W = F,(A) io-o-4^»(^-^) '^'W, (1) 
Es{B -V) = {O.U±0.03)Eg{B -V), (2) 



k'{X) = 2.659 (-1.857 + 1.040/A) + 4.05 0.63 jjm < X < 2.20 /im 

= 2.659 (-2.156 + 1.509/A - 0.198/A2 + 0.011/A^) + 4.05 

0.12 i2m<X< 0.63 fim. (3) 



(4) 



Fo(A) and Fj(A) are the dust-obscured and intrinsic stellar continuum flux 
densities, respectively; A' (A) is the dust obscuration suffered by the stellar 
continuum; E5(B— V) and Eg(B— V) are the color excess of the stellar contin- 
uum and of the nebular emission lines. The latter is measured using hydro- 
gen line ratios and foreground dust ( [Calzetti, Kinney fc Storchi-Bergmann 



1996| ). k'{X) is the starburst obscuration curve, and (3 is the UV continuum 
slope measured between 0.125 fim and 0.26 fim. /5 is a sensitive indicator of 
both reddening ( Palzetti, Kinney fc Storchi-Bergmann, 19"94D and obscuration 
( |Meurer, Heckman fc Calzetti, 1999| ; [Calzetti et al., 200U| ), owing to the fact 



the the UV spectral shape is fairly constant when young, massive stars are 
present. The value of the intrinsic UV slope, Po, depends slightly on assump- 
tions on the star formation history of the starburst; /5o=— 2.1 and —2.3 for 
constant star formation over 1 Gyr and over 100 Myr, respectively. 

Equations 1-4 are valid for all starbursts with optical depths less than a 
few, typically E(B— V)^ <1. In the specific case of Figure 3, we measure 
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Fig. 3. The UV-to-FIR SEDs of two low redshift starburst galaxies: the 
FIR-luminous NGC6090 (left panel) and the Blue Compact Galaxy Toll924-416 
(right panel). UV-optical spectrophotometry is reported as a continuous line in the 
wavelength range 0.12-1.0 ^m, nearlR photometry as empty circles, IRAS data as 
filled squares, and ISO measurements as filled triangles ( Calzetti et al., 200(]| ). Up- 
per limits have downward arrows. The smooth curves represent the best fit to the 
data longward of 40 /im, with a dust emissivity oc and multiple temperature 
components. 

E(B-V)g=0.60 in NGC6090 and E(B-V)g=0.02 in Toll924-416. The for- 
mulae above can be applied to the SEDs of starbursts to calculate the amount 
of ionizing and non-ionizing stellar energy 'lost' to dust absorption, and there- 
fore predict the FIR dust emission. A comparison with the observed FIR 
emission in Local starbursts shows that the agreement between predictions 
and observations is within a factor of 2 for each individual object, and within 
~20% on a sample ( [Calzetti et al., 2000|) . This implies that Equations 1-3 
effectively recover the intrinsic stellar light. The ratio of the observed fluxes 
in the 1-1000 fim range, F(FIR), and in the ultraviolet at 0.16 /xm, F(UV), 
is a measure of the total dust obscuration suffered by the starburst, and can 
be expressed in terms of the other dust-obscuration indicator, /?, as ( [Vleurer 
Heckman fc Calzetti, 1999| ; [Calzetti et al., 2000|) : 



F{FIR)/F{UV) = 1.58[10(°-^2''^+^-9^) - 1] 



(5) 



The wavelength coverage of the FIR SEDs used in this analysis is much broader 
than in many previous studies, thanks to the addition of ISO data; the data 
span the range ~8-240 /im (IRAS-I-ISO), sampling well beyond the peak of 
the FIR emission and probing the contribution of dust cooler than ~30 K. It is 
useful to compare the flux detected in this broader range with the 40-120 /xm 
flux, which has been used to characterize the FIR emission from galaxies when 
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only IRAS data were available: 



F(8 - 240) 

^ ' - 1.72 ±0.25. (6) 



F(40 - 120) 



For starbursts, the total FIR emission is almost a factor of 2 larger than the 
flux detected in the 60 /xm and the 100 /im IRAS bands, and the flux beyond 
120 /im represents ~20-25% of the total QCalzetti et al., 200q ). 



In addition to providing a simple and purely empirical recipe for dust obscura- 
tion corrections, Equations 1-5 give us insights into the physical conditions of 
UV-selected starbursts. Equations 1 and 5 imply that the dust obscuring the 
starburst behaves as if it is external to (surrounding) the region(s) of active 
star formation. A starburst environment is likely to be inhospitable to dust, 
due to the high energy densities, to the presence of shocks from supernovae 
explosions, and to the possibihty of gas outflows. Equation 4 means that the 
stars ionizing the gas are statistically well mixed with the stars producing the 
UV continuum, although Equation 2 implies that the ionized gas suffers about 
twice as much reddening as the stars. A way out from this empasse is to as- 
sume that the ionized gas is more closely associated to the dust than the stars 
responsible for the U V emission (|Calzetti, Kinney fc Storchi-Bergmann, 19941) . 



The ionizing stars may or may not be as dusty as the ionized gas ( |Calzetti" 



1997| ; [Chariot fc Fall, 2000| ), while the non-ionizing, UV-emitting stars will gen- 
erally be less obscured than the gas ( |Ualzetti, Kinney fc Storchi-Bergmannj 
1994| ), as they live long enough to diffuse into regions of lower dust density 
( Garmany, Conti, fc (Jhiosi, 198^ ). HST imaging of nearby starbursts shows 



that both scenarios are present, i.e., ionized gas and dust can be both phys- 
ically separated from and associated with the exciting stars ( |lVlaiz-Apellaniz| 



fc Walborn, 200"0| ). In either scenario, the nebular emission from the galaxy 
will be more absorbed than the spatially-integrated stellar continuum, hence 
Equation 2. 

Equations 1-5 are generally not applicable to the most FIR-luminous objects, 
i.e. ULIRGs. In these galaxies, the gas and dust involved in the starburst are 
centrally concentrated with densities ~10 times higher than those of Giant 
Molecular Clouds (|Downes fc Solomon, 1998)) . Dust is heavily mixed with 



the newly formed stars, and the dominant geometry is very likely to be more 
extreme than in UV-selected galaxies. Recent results from HST show that the 
application of Equation 5 to ULIRGs underestimates the FIR flux by factors 
RilO-100 ( IGoldader et al., 200q ). 
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4 Measuring Dust in UV-Selected High Redshift Galaxies 



The population of Lyman-break galaxies at redshift z>2 ( pteidel et al., 199^ , 
1999| ) presents a number of observational properties which resemble the central 



regions of Local UV-selected starburst galaxies. The high-redshift objects have 
typical SFRs per unit area ?»l-2 yr~^ kpc~^ ( [Meurer et al., 19^ |Calzett 



fc Heckman, 199"9| ); the restframe UV spectra show the wealth of absorption 
lines which is generally observed when young, massive stars dominate the 
emission; nebular emission is also observed in the restframe optical spectra 
(pettini et al., 1998| ), confirming that the Lyman-break galaxies are actively 



forming new stars (i.e., they are not in a post-starburst phase). Finally, the 
range of UV slopes of the high-redshift galaxies is reminiscent of the one 
observed in the Local, UV-selected starbursts, with values — 3 < /3 < 0.4. 

In Local starbursts the large range of UV slopes is a signature of dust obscu- 
ration, and the same interpretation has been suggested for the high redshift 
galaxies, based on the similarities presented above. The median UV slope of 
the Lyman-break galaxies is /3 ~ —1.4, corresponding to about 80% of the UV 
light at 0.17 fim being reprocessed by dust in the FIR. Incidentally, /5 ~ —1.4 
is also the median value of the Local, UV-selected galaxies; the analogy should 
not be taken any further as different selection biases are at work in the two 
samples. The peak of the FIR emission from the z^3 galaxies should be ob- 
servable at wavelengths ^^^200-400 /im (restframe ~50-100 fim). Presently, the 
most sensitive sub- mm instrument, SCUBA at the JCMT, reaches the faintest 
limits at 850 fim, i.e. restframe wavelength f«200 /xm. This is well beyond the 
peak FIR emission observed in galaxies in the Local Universe (e.g.. Figure 3). 
The FIR energy detected at such long wavelength represents a fairly small 
fraction of the total dust emission; in the specific case of the two starbursts in 
Figure 3, the energy emitted at 200 fim represents between 12% (NGC6090) 
and 1.3% (Toll924— 416) of the total FIR energy. Reconstructing the dust 
bolometric emission of high redshift galaxies from SCUBA detections is thus 
bound to suffer from large uncertainties, up to a factor ^10. 

If the two galaxies in Figure 3 were placed at z=3, the fluxes observed at 
850 would be 0.92 mJy and <0.01 mJy, respectively, for NGC6090 and 
Toll924— 416. The R band magnitudes would be 25.7 and 26.4, respectively, 
about 1.5 mag fainter than the ~R* value determined for the Lyman-break 
galaxies (|Steidel et al., 199"9| ) . The two galaxies differ by less than a factor of 2 



in the restframe UV luminosity, and by about a factor of 400 in the restframe 
200 fim luminosity. In addition to have a FIR emission which is about 50 
times fainter, Toll924— 416 has also a much warmer dust SED than NGC6090 
( |Calzetti et al., 2000| ). This explains why in the long- wavelength tail of the 



dust emission the discrepancy in the intrinsic flux between the two galaxies 
exharcerbates. 
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If the luminosity of a z=3 galaxy with the SED of NGC6090 is boosted up to 
be detectable by SCUBA (adopting as fiducial value the 5 a limit of 2 mJy in 
the Hubble Deep Field, [Hughes et al. (199"8| )), the galaxy would have R~24.9 
and a bolometric luminosity Lboi=l-5x 10^^ L©. Thus, such a galaxy would 
need to be as luminous as a ULIRG to be observed with SCUBA ( [Lilly et aT 



1999| ). A more extreme situation is seen with Toll924— 416, due to the very 
small amount of energy emitted at 200 yum, less than 1% of the bolometric 
luminosity. To be detectable with SCUBA, a z=3 galaxy with the SED of 
Toll924— 416 would need to have a bolometric luminosity \j\joi ~3xl0^^ Lq, a 
few times brighter than the brightest non-AGN powered galaxy known in the 
Local Universe. Such a galaxy would have R~19, almost 4 mag brighter than 
the brightest Lyman-break galaxy identified so far (Giavalisco, 2000, private 
communication). Clearly, a galaxy with a luminosity within observed limits 
and with a FIR SED as warm as that of Toll924— 416 would not be detectable 
with SCUBA. 

To summarize, a z=3 galaxy with an SED similar to that of Local UV-selected 
starbursts needs to be dusty, as bright as a ULIRG, and to have a relatively 
cool FIR SED to be detectable with current sub-mm instrumentation. The 
first two criteria already restrict the number of candidates, since, even if high 
redshift galaxies may be dusty, it is not established how many of those have 
bolometric luminosities in the 10^^ L0 range. The third criterion (cool FIR 
SED) is equivalent to require that, for an emissivity index e=2, dust with 
temperature ~20 K must be present in the galaxy; this criterion may further 
bring down the number of candidates, depending on the details of the distri- 
bution of dust and star formation within the galaxy. A recent project aimed 
at observing Lyman-break galaxies with SCUBA has yielded only one ten- 
tative detection out of ~10 targets ( [Chapman et al., 2000|) , confirming that 



a number of factors have to play together to make the long wavelength tail 
of the FIR emission of those galaxies detectable. Photometric uncertainties 
in the restframe UV measurements and Malmquist bias add to the difficulty 
of making reliable predictions for the FIR luminosity of the high-z galaxies 
([Adelberger fc Steidel, 2000| ). 



The effects of dust in the distant Universe still need to be placed on firmer ob- 
servational ground. To achieve this target, high redshift star-forming galaxies 
require better spectral energy coverage. The available optical observations will 
need to be supplemented with near- to mid-IR spectral data to measure the 
restframe optical emission and nebular lines, and with far-IR/sub-mm/mm 
data to derive a full picture of the dust emission. SIRTF, ALMA, and NGST 
are the key missions which will address this issue. Knowing how much dust is 
present in galaxies at different redshifts will not only allow us to derive the 
intrinsic star formation history of the Universe, but will also shed light on a 
number of other astrophysical issues, such as the dust and metal enrichment of 
the evolving galaxies and of the Universe, and the starburst-AGN connection. 
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